Phase separation transition in anti-ferromagnetically interacting particle systems 
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One dimensional non-equilibrium systems with short-range interaction can undergo phase tran- 
sitions from homogeneous states to phase separated states as interaction (e) among particles is 
increased. One of the model systems where such transition has been observed is the extended Katz- 
Lebowitz-Spohn (KLS) model with ferro-magnetically interacting particles at e = 4/5. Here, the sys- 
tem remains homogeneous for small interaction strength (e < 4/5), and for anti-feromagnetic inter- 
actions (e < 0). We show that the phase separation transitions can also occur in anti- ferromagnetic 
systems if interaction among particles depends explicitly on the size of the block (n) they belong 
to. We study this transition in detail for a specific case e = S/n, where phase separation occurs for 
8 < -1. 

PACS numbers: xxx xxx 



One dimensional driven diffusive systems 0, Q have 
been a subject of extensive studies in recent years. Some 
driven diffusive models with local dynamics exhibit ex- 
otic phenomena like phase separation and phase transi- 
tions [3 8] , with novel spatial correlations in their steady 
states, which can not be seen for systems in thermal 
equilibrium, with short range interactions. Since there is 
no general theoretical method to study non-equilibrium 
systems, phase transitions for them can not be inferred 
from any guiding principle. Few years ago a general cri- 
terion for the existence of phase separation in density- 
conserving one dimensional driven diffusive systems, have 
been proposed by using a mapping of these models to a 
Zero Range Process (ZRP) by Kafri et al. [|. It has 
been suggested there, that existence of phase separa- 
tion in any given model depends only on the rates (or 
the steady state curret J n ) at which domains of vari- 
ous sizes (n) exchange particles. A strong phase separa- 
tion, where density fluctuations are limited to the domain 
boundaries, occurs for arbitrarily small densities when 
the steady state current J n flowing through a block of size 
n vanishes for thermodynamically large blocks. Whereas 
a condensed phase with large density fluctuations occurs 
when the current 

J„~ Joo(l + b/n a ) (1) 

with either a < 1 and b > or with a = 1 and b > 2. In 
this case phase separation occurs only for large enough 
densities. 

This criteria [9( is based on a mapping which relates 
driven diffusive systems with a well known model called 
Zero Range Process. In ZRP, each site on a lattice (for- 
mally called boxes) can accommodate more than one par- 
ticle, which are allowed to hop to the neighboring box 
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FIG. 1: A microscopic configuration of driven diffusive system 
is mapped to corresponding configuration in ZRP. 



with a rate u(n) that depends on the number of par- 
ticles n present in the departure box (see Fig [1} . The 
vacant sites of any particular driven diffusive system can 
be considered as the boxes and the domain of size n to 
the immediate left of this vacant site can be considered as 
the number of particles. Finally, the steady state current 
J n can be regarded as the hop rate uin) of ZRP. Such a 
mapping can be shown to be exact [9( in some models. 
In others [f| 0], though the ZRP mapping is only ap- 
proximate, phase separation could be predicted beyond 
doubts. Phase separation transition with a = 1, b > 2 
was first observed in a class of models Q where two 
species of hard-core particles (equal in number) on a one- 
dimensional ring interact ferro-magnetically, which was 
extended later [7j for unequal density of particles. In 
these cases, b depends on the interaction strength e and 
condensation occurs for large density of particles as e is 
increased beyond a critical value e c . For e < e c , includ- 
ing the antiferromagnetic interaction e < 0, the system 
remains homogeneous. 

In this article we introduce a model with two species of 
particles which interact antiferromagnetically, where the 
interaction strength depends on the size of the domain. 
We show that such modifications can lead to phase sep- 
aration transition, even for anti-ferromagnetically inter- 
acting particles. An intuitive and physical argument is 
provided in favour of the possibility of such a transition. 

The model is defined on a ring where sites are la- 
beled by i = 1,2, ...L. Each site i can either be vacant 
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(si = 0) or occupied by a positive (sj = +1) or a neg- 
ative (si — —1) particle. These particles interact anti- 
ferromagnetically with interaction strength e(nj) which 
depends on the size ni of the I th block that the particles 
belong. Note that a block (or a domain) here is defined by 
an uninterrupted sequence of particles and leveled by an 
index I = 1, 2 . . . No, where No is the number of vacancies 
in the ring. Formally the interaction can be represented 
by 
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The model evolves according to the nearest-neighbor ex- 
change rates 
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-0. (3) 



Clearly, the dynamics is particle conserving. Thus the 
system can be characterized by the relative densities p = 
^±±^ an d „ = 



where N± are the number of 



N N + +N- 

and — particles respectively. 

First let us consider the non- interacting case, e = 0. 
Here, the steady state weight 14 7 ({&;}) of all the con- 
figurations in grand canonical ensemble (GCE) can be 
written as Is 01 
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W({m}) = l[z^z ni , 



(4) 
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where ni is the number of particles that reside in the 
I th block and Z ni is the sum of all microscopic config- 
urations of a block of size n.;. The fugacity z is associ- 
ated with both (±) kinds of particles. Since the weight 
of the configurations are factorized in terms of weight 
z n 'Z ni of individual blocks, one can draw an analogy of 
the model with a ZRP having No boxes and ( N + + N- ) 
particles where steady state in GCE has a product mea- 
sure; the single box weight of ZRP is f(n) = Z n . Fur- 
ther, from the exact results of ZRP in grand canonical 
ensemble f(n) = J\k=i one § ets u ( n ) = Zn-i/Zn 

(as Z n — rifc=i Zk/Zk-i)- Since in the non interact- 
ing system, the current through a block of size n is 
J n = Z n -i/Z n , one identifies u(n) = J n . Now, from 
the assymptotic behaviour § of J n = \ (1 + ^ ) one can 
conclude, using criteria mentioned in Eq. (fT]), that phase 
separation is not possible in this case. 

The exact mapping of the model to ZRP does not hold 
for interacting system. However, it has been argued @ 
that the correlations between neighbouring boxes can 
still be neglected for sufficiently large a and the corre- 
spondence u(n) = J n provides a correct prediction of 
phase separtion transition. Since the blocks are consid- 
ered uncorrelated, one can calculate J n by modeling a 
single block of size n where particles follow dynamics ([3]) . 
Steady state weights of such a system on a ring [H has 
an Ising measure P({si — ±}) exp(—/3H) with 
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i.e., weight of any configuration is same as that of the 
above system in equilibrium. The second term is needed 
to fix magnetization of the Ising system M. — 2j]—l. For 
equal density case, h — 0. Thus, the current 
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is a static correlation function of the Ising system Eq. (|5|) , 
which can be calculate d II 01) using the standard trans- 
fer matrix formalism [111 ]. For a thermodynamically 
large block, current Joo calculated on a ring is identi- 
cal to that of an open domain. However, the corrections 
Jn = JoaO- + b(e)/n) depends on the boundary condi- 
tions. It has been argued [l2j that 6(e) for an open sys- 
tem is related to that of the ring by a universal factor; 
6(e) = where the subscript R stands for ring. Thus 
it is sufficient to calculate correlation functions Eq. © 
on a ring with respect to the equilibrium Ising measure. 
A simple transfer matrix formalism reveals, 
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FIG. 2: Phase diagram in w-S plane: the shaded region cor- 
responds to phase separated state. Inset, here, shows depen- 
dence of b on e for equal density case. 

These relations predict that, for any relative density 77 
the system undergoes a phase transition from a homoge- 
neous state to a phase separated state when e is increased 
beyond a critical threshold e c . Inset of Fig. [2] shows vari- 
ation of 6(e) with e for equal density case i] — 1/2; here 
e c = 4/5. It is evident from the figure, that phase sep- 
aration is not possible when e < (anti-ferromagnetic 
interaction) . These predictions for ferro-magnetically in- 
teracting system have been verified using Monte Carlo 
simulations of extented KLS (EKLS) model for both 
7? = 1/2 @ and for 77 ^ 1/2 0. 

A possible reason, why antiferromagnetically interact- 
ing particles do not phase separate is the following. For 
ferro-magnetic case, each particle likes the same kind 
of particles as neighbours and the asymmetric exchange 
H — — > — h, then arrange more + particles towards the 
right end of the block. With this arrangement of par- 
ticles, 0s enter easily into the ferromagnetic block from 



both sides allowing the movement of the whole block as a 
single entity. However, in an antifcrromagnctic system, + 
and — particles prefer to be arranged alternatively which 
restrict Os from the left (right) to move further once it 
encounters a + (— ) particle. The antiferromagnetic inter- 
action could compete with the boundary dynamics when 
e ~ 0(l/n), [i.e. ne ~ 0(a)\. Thus, it is suggestive 
that a phase separation could be possible when e„ = — , 
explicitly depends on n. In the rest of the article we 
concentrate on antiferromagnetic case with interaction 
strength e n = — , that depends explicitly on the size of 
the block. It is not difficult to verify that the steady 
state have Ising measure. Thus, J n can be calculated in 
a straightforward way following the procedure discussed 
here. This results in, 

j n = Joo (1 + 3 ~2n J °° ) ; J - = ^ - v)- (8) 

Thus, b(5,r)) = (3 — 45 Joo) /2 and for any given 77 phase 
separation occurs for systems having large density p when 
5 < —[477(1 — t?)]" 1 . Corresponding phase diagram in t]-5 
plane is shown in Fig. [2J 
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FIG. 3: Comparison of density profiles: solid line corresponds 
to the density profile of a block of size n — 100 and L = 500. 
Symbols correspond to that of a ring with 100 particles and 
a single defect 0. In both cases, S = —0.9. 

To check the correctness of the prediction, we choose to 
study only the equal density case in details and verify the 
predictions there. In this case, = 1/4 and b(S) = 
thus the transition is expected for S < — 1. Note, that 
in deriving Eq. ([5]) we have assumed that, a domain con- 
taining n particles of any species, behaves like an open 
system where (+) particles enter from left and (— ) par- 
ticles exit from right with rate a. Hence, the dynamics 
of a domain can be modelled by a ring having n particles 
(±) and a single defect 0, and following above dynamics. 
To check the validity of this assumption we have calcu- 
lated the density profile of a block of size n = 100 using 
Monte Carlo simulation of a system of size L = 500 at 
S = —0.9. In Fig. [3] the density profiles of a block of 
size n = 100 obtained from Monte Carlo simulation of 
the full system is compared with that obtained from a 
ring with 100 particles and a defect 0. In both cases we 
choose 5 = —0.9. An excellent match seen here, justifies 



the assumption that a domain of size n in fact behaves 
like a ring having n particles and a single defect 0. 




FIG. 4: Block size distribution for S = -0.2 (black) and -2.0 
(red online) for a system with L — 2000 and No = 100. 

The coarsening dynamics of EKLS model is quite slow. 
Thus in fact achieving a true condensate through Monte 
Carlo simulations, in a finite time, is practically impos- 
sible. However, the distribution of block sizes P(n), pro- 
vides an indication: when 5 < 5 C = — 1, a condensate 
appears along with the usual scale-free distribution. In 
fact, to observe the condensate one must take the den- 
sity much larger than p c = 1/(6(5) — 2). Figure @] shows 
distribution of size of the block for 5 = —0.2 and —2.0, 
correspondingly the values of 6 are 6 = 1.6 and 2.5. The 
relative density is given by p — 19 with L = 2000. It is 
clear that condensation emerges only in the latter case. 
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FIG. 5: Evolution of a phase separated configuration: the 
snapshots of the system are taken at t = 0, 10 3 , 10 4 , and 10 6 . 
(a) Upper and (b) lower pannels correspond to S = —0.2 and 
-2.0 respectively. Here L = 2000 and iV = 100. The double 
arrow (red) corresponds to a large condensate that survives. 

It is expected that, in the phase-separated regime, any 
initial phase-separated configuration evolves only to gen- 
erate fluctuation in the size of the condensate; the con- 
densate does not disappear. To check the stability of 
an initially phase separated configuration we have done 
Monte Carlo simulation of the system of size L = 2000 
with No = 100 and look at the space-time plots for 
2000 MCS at t = 0, 10 3 , 10 4 and 10 6 . The snapshots 
for S = —0.2 and —2.0 arc shown in Fig. \5\ Clearly for 
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the latter large condensate (indicated by a double 

arrow) survives even at t = 10 6 . 

In conclusion, we have studied the extended KLS 
model with antiferromagnetically interacting particles. 
The strength of interaction between particles e(n) — 5/n 
in this model depends explicitly on the size of the block n 
they belong to. The model can be mapped to zero range 
process by identifying vacancies as boxes; the hop rate 
u(n) in ZRP is identified as the particle current J n in the 
lattice, which can be calculated exactly. Such a mapping 
enables one to predict that the phase separation transi- 
tion (same as condensation transition in ZRP) occurs for 



5 < — 1. These predictions are verified by using Monte 
Carlo simulations. 

Some comments are in order. Unlike the particle cur- 
rent in the lattice, the hop rate in ZRP is stochastic and 
uncorrelated in time. Recent studies [l3| of ZRP with 
time-correlated hop rate suggest that effective b which 
sets the criteria for condensation transition may change 
when the correlation is significantly large. However, our 
numerical study of the antifcrromagnetic system indi- 
cates that the time-correlations, present in J n , are pos- 
sibly irrelevant. In fact, such correlations are known Q 
to be irrelevant in AHR [4j model. 
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